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ABSTRACT 

Fireball model of the GRBs predicts generation of numerous internal shocks, which 
then efficiently accelerate charged particles and generate magnetic and electric fields. 
These fields are produced in the form of relatively small-scale stochastic ensembles 
of waves, thus, the accelerated particles diffuse in space due to interaction with the 
random waves and so emit so called Diffusive Synchrotron Radiation (DSR) in contrast 
to standard synchrotron radiation they would produce in a large-scale regular magnetic 
fields. In this paper we present first results of comprehensive modeling of the GRB 
spectral parameters within the fireball/internal shock concept. We have found that 
the non-perturbative DSR emission mechanism in a strong random magnetic field is 
consistent with observed distributions of the Band parameters and also with cross- 
correlations between them; this analysis allowed to restrict GRB physical parameters 
from the requirement of consistency between the model and observed distributions. 

Key words: acceleration of particles - shock waves - turbulence - galaxies: jets - 
radiation mechanisms: non-thermal - magnetic fields 



1 INTRODUCTION 

The fireball model is currently accepted as a standard model of the gamma-ray burst (GRB) prompt emission (e.g.,[Meszaxog 
20061 ) . It is supposed that a central engine produces a number of relativistic internal shocks, which then interact with each 
other. The phenomenon of the shock waves requires an efficient mechanism of energy dissipation. In a coUisionless case, the 
most efficient ways of the energy dissipation are via generation of fluctuating electromagnetic flelds and acceleration of charged 
particles up to high energies. 

Microscopically, this field generation can be d r iven by two-stream instabilities associated wi t h the shock propaga- 
tion (iKazimura et alJ[l998l: iMedvedev fc Loeblll999l: [JYederiksen etlo] |2004|: Isret et aLllJooi. boosi: JNishikawa et aflbood: 



Jaroschek et al.ll2004l : lHededaj2005l : lHededal fc Nishikawal2005l : lBret fc Dieckmannl2008l : lKeshet et al.l2008l : lDieckmann fc Bret 
2010|), while the ac celeration of pa r ticles is provided by their interaction with the shock-gene rated random and regular elec 



tromagnetic fields (JMeszarosI [2OO2I : iNishikawa et aflboosl : IPiranlboOSi : ISaril bood : [silvalbood ). It is well established by now 
that the magnetic and electric fields produced in the shock interactions have often a significant random component at various 
spatial scales. 

The presence of the random component is critically important for generation of nonthermal emission from corresponding 
objects. Indeed, unlike regular gyration in the presence of a regular magnetic field, the shock-accelerated charged parti- 
cles moving through a plasma with random electromagnetic fields experience random Lorenz forces and so follow random 
trajectories representing a kind of spatial diffusion. Accordingly, the particles produce a diffusive radiation whose spectra 
depend on the type of the fie l d (magnetic or electric) and on spectral energy distribution of t he field over the spatial scales 
(iToptygin fc Flei shman 'l987l: lilededall boosi : IPleishmanl bood : [Fleishman fc BietenholjIioO?! : [Fleishman fc ToptygiiilboOTEil : 
Sironi fc Spitkovskv..2009. V Be low in this paper we rely on analytical DSR theor y proposed bv lToptygin fc FleishmanI ( 19871 ) 
and then further developed bv lFleishman| ( | 2006|) : Fleishman fc Bietenhola ( 2007 ). An alternative way of calculating radiation 
is the use of numerical PIC simulations l)Hededalll2005l : lHededal fc NishikawalbOOa : INishikawa et al.ll2008l : ISironi fc Spitkovskv 
which confirm the analytical results in the common parameter domain. However, the case of strong random field and 
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strong angular scattering of the radiating electrons requiring a large dynamic range of the involved parameters is yet beyond 
available PIC capacities, which justifies the choice in favor of the well tested analytical theory. 
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Individual spectra of the prompt GRB emission are typically well fitted by a phenomenological Band function IjBand et al. 



19931 ) ■ which consists of low-energy (spectral index a) and high-energy (spect ral index /3) pow er-law regions smoothly linked 
at a break energy Etr- The diffusive synchrotro n radiation (DSR ) was shown ( Fleishman 20061) to produce spectra consistent 



with that observed typically from the GRBs IJBand et al. 



I993I : 



Mazets et alJ 



2004 



Ohno et al.l boosi : IPal'shin et aD boOSi : 



Granot et al.ll2009l ). It is y et unclear, howeve r ^ if the DSR spectra are naturally consistent with observed distribution of the 
GRB spectral parameters ( Preece et al.ll200d : lKaneko et al.ll200a ) and what ranges of physical GRB parameters are needed 
to reconcile the theoretical spectra with the observed ones. In this paper we present a model of GRB prompt emission 
generation by DSR in relativistically expanding GRB jets. The input parameters of the model are constrained by available 
observations and take into account dependences between involved parameters implied by physical laws. We vary a number of 
free parameters of the model to achieve the best agreement between the variety of the modeled and observed spectra. This 
analysis confirms that the DSR model, specifically — the non-perturbative strong-field regime, is intrinsically consistent with 
the observed distributions of the GRB spectral parameters. 



2 FORMULATION OF THE MODEL 

To be specific we adopt the fireball model in which the GRB prompt emission is generated in a coUimated jet ejected with a 
relativistically high speed 11 from a central engine. Adopting a general internal shocks/fireball concept we accept that a single 
binary collision of relativistic internal shocks results in a single episode of the GRB prompt emission. Microscopically, this 
shock-shock interaction first produces high levels of random magnetic and/or electric fields and accelerates the charged particles 
up to large ultrarelativistic energies; and then, these particles interact with the random fields to generate the gamma-rays. 
Although there are some common general properties of all cases of relativistic shock interactions, each shock-shock collision is, 
nevertheless, unique in terms of combination of the physical parameters involved. Accordingly, we are going to estimate and 
adopt a set of standard ("mean") parameters appropriate to account for the most global GRB properties, and then consider 
if a reasonable scatter of those standard parameters is capable of reproducing more detailed properties of the considered class 
of events as a whole — the statistical distributions of the GRB spectral parameters and cross-correlations between them. To do 
so, we consider a number of different emission models including the standard synchrotron radiation and DSR regimes in case 
of either weak or strong random magnetic field. The spectral slopes and breaks depend on both the emission mechanism and 
combination of physical parameters affecting the radiation spectra within a given mechanism. Thus, the goal of the modeling 
is to establish if there exists a parametric space making one or another theoretical model compatible with the observational 
data on the GRB spectral properties. 



2.1 Basic parameters of the GRB source 

From analysis of so-called 'compactness problem', it is established that the bulk Lorenz -factor of the expanding jet F 



1/^/ 1 — v'-^/c^, where c is the speed of light, must be much larger than unity, F > 100 ( Paczvnski fc Rhoadd 119931 : IPiran 



2OO5I : ISarilbOOq ) . A maximum F value is not well constrained observati onally; being co nservative, we will not consider values 
above 1000. The total kinetic energy of the jet is roughly E ~ 10^^ erg I Meszar os 2006), which, along with F estimate and an 
assumption of particle composition, allows estimating of the total number of ejected particles. We adopt that bulk of the jet 
mass resides in the protons, thus, the total number of ejected protons is estimated as 

Af = S/(FmpC^) ~ 10^\ (1) 

where rUp is the proton (rest) mass. If no e+e~ pairs are produced, then the number of electrons is equal to the evaluated 
number of protons. 

The number density of the particles is 

n=^, (2) 

where V is the jet volume. The full volume occupied by the jet material is a cone with the opening angle Q,ang (for the 
modeling we adopted a constant value of Jiang = 10~^ ster), so Vcone = (47r/3)r2angi?'^, while the volume of spherical layer 
with thickness A of this cone is 1/ = {ATv/3)Qang{R'^ — (R ~ A)'^). The volume is not well constrained by the observations, 
although it can be estim ated using observed time scale St of the emission variability. Indeed, for a relativistically moving 
source (e.g.. |Piranll2005l ) we have St ~ A/(2cF^); accordingly, for the co-moving frame we can estimate: 
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Then, the number density in the co-moving frame is 



(3) 



(4) 
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Figure 1. Model histogram of the spectral index distribution assuming slow cooling synchrotron regime. Note, that most of the indices 
are displaced by a small value compared with the asymptotic value, —2/3 (the 'line of death'). There are a number of outliers with a few 
of them apparently violating the 'line of death'. These outliers originate from the fit errors provided that the Band fitting function does 
not represent a perfect match to the theoretical synchrotron spectrum. 



The radiation spectra for each set of the source parameters are calculated in the co-moving system and then the emission 
frequency is transformed to the observer's frame taking into account the relativistic motion of the source and cosmological 
expansion of the Universe 



^'^^(1 + .). 



(5) 
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Let us introduce the energy contents (e.g.. lSarill2006f ) of constituents needed to produce radiation — the magnetic field 
Ub B^ 



8tvw 



where w = Vnmp(? is the kinetic energy density of the expanding shell and Ub 
field; and accelerated electrons 

Ue -y^errie 
w 1 rrir, 



(6) 



is the energy density of the magnetic 



(7) 



where Ue = Ue'y^emsC? is the energy density of accelerated electrons, f e is a fraction of electrons being accelerated (^e < 1 if 
only a fraction of all electrons are accelerated ( Bvkov fc Meszarodll99a ). while others remain "cold", and ^e > 1 if e^e~ pairs 
are produced at the emission source), 7 is a characteristic Lorenz-factor of the accelerated electrons, and me is the mass of 
electron. 



2.2 Standard synchrotron model 

Although synchrotron models are generally consistent with overall G RB energetics and lig ht curves, they are intrinsically 
incompatible with the distribution of low-energy spectral index a (e.g.. lBaring fc Brabvll2004l ). To demonstrate this explicitly, 
we show the model distribution of the low-energy indices obtained within a synchrotron model in Figure [l] 

To be specific in generating this histogram we assumed the slow cooling regime and reasonable statistical distributions 
of the relevant parameters (see below for greater detail), such as bulk Lorentz factor of the relativistically expanding shell, 
energy content of the accelerated particles and produced magnetic field, number density of the particles, and their e nergy 



spectrum, as well as correlations between these parameters consistent with observations ( Meszaroa 120061 : ISaril 120061 ) . The 
histogram in Fig. [T] rep resenting an asymmetric narrow distribut i on peaking around a = —0.7, is in evident contradiction 



with the observed one ( Preece et al.l l2000l : iBaring fc Brabvl l2004l : iKaneko et al.l 120061 ). which is a more or less symmetric 
broad distribution peaking at a ~ — 1. The fast cooling regime results in a histogram similar to that in Figure [l] with the only 
difference that it peaks around a ~ —1.5. We conclude that a more sophisticated modeling is needed to achieve reasonable 
agreement between the observations and the theory of electromagnetic emission in the GRB sources. To address this problem, 
Medvedev (2000) proposed that emission of fast electrons moving in small-scale random magnetic field may possess the 
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spectral properties consistent w ith those observed from GRBs; the corresponding DSR process in the GRB context has than 
been studied quantitat i vely b y iFleishmanI (|2006l ) within general concept of the stochastic theory of radiation proposed by 
Toptygin fc FleishmanI l| 19871 ). 



2.3 Theory of DSR: main equations and parameters 

For the purpose of this more detailed modeling we note that at the sites where the internal shocks interact in the GRB 
sources, charged particles are accelerated an d two-stream ins tabilities produce high level of random magnetic and/or electric 
fields, so the diffusive synchrotron radiation ( Fleishmanll2006l ) is expected to be produced there. Let us remind basic equations 
describing the DSR and main parameters determining its spectrum. We adopt that the energy of the random magnetic field 
is distributed over spatial scales according a power-law: 



Kai3iqo,q) = -K{qo,ci) [ S^fj 



where 



Kid) 



A, 



A, 



v-l 






[Bi) 



-K{(i)S{qo - <7o(q)) iSc^fi - -^ 



('-~l)Ch.'(ga) 



Atv 



Q'min ^ Q ^ Qniax , 



u is the spectral index of the random field distribution, the spectrum K{q) is normalized by d^q so that 

/<(q)d' q = (-Bst), gmin < gmax, U > 1, 



(8) 



(9) 



(10) 



{B%) is the mean square of the random magnetic field. 

For the adopted spectrum of the random magnetic field, the perturbative DSR spectrum I Fleishmanll2006l ) can be obtained 
analytically: 



Sttc 
where 7^ = 7 



■q{uj) 



(11) 



2 / 



q(uj) = 



i>(i,+2Y^(v+3) 



for Lo < 2^ or oj > 27^(^0, 



7r ai/ijj^ 






:;^ + 






[2^c)\ fo^' 2SI 



< OJ < 27 (jJo. 



(12) 



2vuiQi'-' [^ "^ 4(u+2p"/iui'^ 2(^+3)76 cjg ' 4(i/+2)7*t^g '" \^27jfeocy 

These equations describe the DSR from particles moving along almost rectilinear trajectories and so valid only for 
relatively weak magnetic field. In the internal shock interactions, how ever, a strong random magnetic field can often be 



produced, which requires a more general, non-perturbative treatment ( FleishmanI 120061 : iFleishman fc Bietenhola 120071 ) 
which the DSR spectrum has the form 



lu: = 



q{uj)^2 



Sttc 



2 2 

^pc7 



^(s) 



(13) 



where $(s) = 24s^ /q°° dtexp(-2st) sin(2st) [cotht- 
^ ( 'it) ) ( 1 ^ — 17^~ ) ' ^^'^ scattering rate q{uj) has been defined by Eq. (fT2)l 



is the Migdal's function, which depends on a single parameter s 
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Therefore, the DSR intensity depends on the following (microscopic) parameters: the electron plasma frequency uipc, the 
(defined by the random field) gyrofrequency Wst, the frequency ujo — gminC specified by the main scale Lmax ~ 27r/gmin of the 
random field spectrum, the spectral index v of the random field, the Lorenz-factor 7 of the emitting relativistic electrons, and 
also on rmmber (and spectral distribution) of the emitting electrons. 



2.4 Links between microscopic and macroscopic parameters of GRBs 

One of the difficulties in direct modeling of the GRB spectra is that they depend on microscopic parameters, which are 
basically unknown. However, we ca n link many of them with t he macroscopic and phenomenological parameters of the GRB 
source introduced in ij 12.11 (sf, e.g.. lKumar &: McMahonll2008l ). 

The electron plasma frequency is specified by the electron number density rie 



Airn^e^ 



(14) 



Here we have to substitute the number density Q found from the estimates of the total number of particles ([T| and the 
source volume ((3|, then 
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(15) 



where e and rrie are the electron charge and mass, nip is the proton mass, 5t is the observed time of the emission variability, 
^ang is the jct Opening angle, and E is the bulk kinetic energy of the shell. Therefore, the plasma frequency is expressed 
through a number of values, which can either be directly observed (like St) or estimated from the observations (like E, T, and 

^ I'ang) • 

Similarly, the stochastic gyrofrequency 

^Bt = „.0 O (16) 



is defined by the random magnetic field value, and so can be expressed via the phenomenological parameter tB describing the 
energy content of the magnetic field: 



'^'' ^ Y mlc^V ^ Y mic2^!^„„9(2c5tr)3 ^ Y -imlc^r^St^ ' ^ ' 

The third parameter of the DSR spectrum having the dimension of frequency is ojo ~ gminC = 27rc/Lmax, which is 
determined by the main scale of the random magnetic field distribution. To parameterize this unconstrained value we introduce 
a dimensionless model parameter A — ujq /uj^t so that 

UJQ — k- LOst- (18) 

Then, A ^ 1 corresponds to a weak random field, while A ^ 1 to a strong random field. 

Finally, the distribution of the magnetic energy over the spatial scales is determined by the spectral index u. We have no 
reliable constraint for this value from the GRB observation, however, using analogy with other astrophysical and laboratory 
cases, we adopt 1 < ;/ < 2. 

Consider now parameters characterizing emitting relativistic electrons. The typical Lorenz-factor 7 of electrons can be 
expressed from Eq. ((Tjl via parameters ee, Co ^-nd T: 

7=^. (19) 

Therefore, a single DSR spectrum is fully determined by specifying the following set of ten involved parameters: E, T, e^, eg, 
A, Ce, St, riang, V, aud z. We turn now to specifying the corresponding parameter ranges. 

2.5 Modeling strategy 

For the sake of further modeling of the GRB spectra within the DSR emission model we note that the spectral parameters 
a, /3, and -Ebreak of the Band fitting function display single-mode distributions. In terms of statistical properties of the GRB 
sources this implies that the GRB physical parameters are taken from the same parent distributions. 

Therefore, we will perform a kind of "global" statistical modeling of the GRB spectra. To do so we produce a large number 
(~5000) of individual DSR spectra, which differ from each other because different combination of input parameters is selected 
for each individual spectrum. Specifically, the ten independent parameters are randomly selected from the corresponding 
parent distributions, while the dependent (derived) parameters are then calculated as described in the previous section. 

The independent parameters, where known, are taken based on observations available and the standard fireball model. 
For example, we adopt the total shell energy E ~ 10^^ erg, the jet opening angle i^ang ~ 10"'^ ster, and F ^ 300 (> 100) to 
have normal distributions around the mean values, while the time variabilit y scale St to have a l og-normal distribution with 



the mean (St) = lO" s and standard deviation a{St) = 10'^'^ s as observed ( Nakar fc Piranll2002l ). Parameters te and es are 
poorly constrained by observations; we kept them much less than unity in all cases, while consider a number of distributions 
for the spectral index i/, which is also unknown parameter of the model. To reduce the number of free model parameters, we 
kept some parameters constant, which have only minor effect on the DSR spectrum shape; namely, we adopted 2 = 2 and 

Ce = l. 

Then, we made a simplifying assumption about the energy spectrum of the shock-accelerated electrons. Although this 
energy distribution is likely to be a broad one (for example, a power-law ne{E) oc E~^), which can easily be taken into 
account, we adopt here a monoenergetic electron distribution. Indeed, adding the power-law energy distribution to the model 
would yield some spectrum regimes common for DSR and other competing emission processes including standard synchrotron 
radiation. We, however, want to evaluate the capability of the DSR mechanism itself to reproduce the Band function parameter 
distributions. If we succeed to do this with the monoenergetic distribution, then the power-law distribution will also be an 
acceptable one, since this will only increase (but not decrease) the variety of the produced radiation spectra; we return to this 
point later. 
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Figure 2. Change in the DSR spectrum due to bulk Lorenz-factor F change. 
Figure 3. Change in the DSR spectrum due to ee change. 

Finally, parameter A is unknown. To evaluate the range of this parameter the most consistent -with observed distributions 
of the GRB spectra, we vary it in a broad limits in our modeling from run to run, while keep a constant within each run. 



3 MODELING RESULTS 

3.1 Dependences of the DSR spectrum on the input parameters 

Before turning to the statistical modeling described, it is worthwhile to consider how the DSR spectra depend on the indepen- 
dent input parameters. These dependences are not self-evident, because each input parameter (e.g., the bulk Lorenz-factor or 
bulk energy of the shell) affect a few microscopic parameters, which, in their turn, specify the DSR spectrum produced. 

To do so, let us select the following set of "standard" parameters F = 300, te = 10~ , e_B = 10~ , A = 2, ^e = 1, (5t = 1 s, 
^ang ~ 10~ , V — 1.2, aud 2 = 2 aud then consider how the DSR spectrum changes as one of these parameters changes, while 
the others are kept the same. 

Figure [2] shows that the DSR level decreases and the spectrum moves towards lower energies as the bulk Lorenz-factor 
increases in contrast to simple expectation based on the bulk doppler shift only. The reason for such a behavior is that the 
bulk Lorenz-factor (for the same other input parameters) affects also other microscopic parameters, so the net effect of the 
F change in the fireball model is different from purely kinematic Doppler effect. Increase of the relativistic electron energy 
content te shifts the whole spectrum towards higher energies. Figure |31 since the main effect of e^ increase is the corresponding 
increase of the typical Lorenz-factor 7 of the relativistic electrons. 

Figure |4] then displays that increase of the magnetic energy density, es, results basically in the upward shift of the whole 
spectrum and also to some shift of the high-frequency part of the spectrum towards higher energies. Parameter A = ojo/ti^st 
has a major effect on the DSR spectrum Fig. [5] because it is the parameter that controls, which DSR regime, perturbative or 
non-perturbative, is in fact realized for a given parameter combination, i.e., weak or strong random magnetic field is present 
at the source. Accordingly, the DSR spectrum shape changes significantly as A changes, new spectrum asymptotes arise as it 
decreases and non-perturbative DSR regime in a strong field develops. 

Decrease of the ^e parameter results in a shift of the spectrum towards higher energies. Fig. [6l The reason for that is the 
same as at te increase: according to Eq. p9[) both of them lead to 7 increase and so to emission of correspondingly enhanced 
(oc 7^) energies. Increase of the cosmological red shift parameter shifts the whole spectrum towards lower energies. Fig. [T] as 
can be expected from the kinematics because no microscopic parameter depends on z in this model. 

Finally, the spectral index v of the random field affects the high-energy slope of the spectrum in the perturbative case, 
while affects other spectrum asymptotes as well in a more general, non-perturbative case, as can be seen from Figure [8] Now, 
when the important dependences of the DSR spectrum on the input parameters are established, we are in the position to 
perform the statistical modeling, and, to fine tune the parameter ranges towards model distributions matching the observed 
ones. 

We start the modeling by adopting trial distributions for the parameters we are going to vary within a single run. 
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Figure 5. Change in the DSR spectrum due to A = uJo/uiBt change. 
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Figure 6. Change in the DSR, spectrum due to ge change. 
Figure 7. Change in the DSR spectrum due to z change. 

Specifically, we initially adopt normal distributions for the involved parameters "with the following means and standard 
deviations: (F) = 300, ar = 50, (te) = 10-^ a,, = 10-^ {^b) = 10-^ a,^ = 10-^ {log^g 5t) = 0, oi^^^^h = 3 ■ lO'S 
(i^) = 1.2, and cr^ = 5 • 10~^. Whenever possible, we will check if other than normal distribution ofi'ers a better fit to the 
observational data. 



3.2 Case I: Week random magnetic field; perturbative treatment applies 

We be gin with the simplest case when the DSR spectrum can be described by a perturbative formula e derived in 



11200611 ■ which are w idely used to evaluate the DSR in the perturbative ("jitter") regime (see, e.g., IWorkman et al 



I200S 



Fleishman 



Mao fc Wand 120071 ). The condition for random magnetic fi eld to be weak is A = cjo/'^st 2> 1; then, applicability of the 
perturbative treatment requires additionally 7 < ajpoCJo/cjgt ( Fleishmanll2006l ) . Adopting a constant value A = 10 and taking 
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Figure 8. Change in the DSR spectrum due to v change. 
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Figure 9. Example of the model noisy DSR spectrum (symbols, A = 2, f = 1.42) and fitting Band spectrum (solid thick green line). Thin 
solid lines display the low-energy and high-energy power-law asymptotes of the fitting Band function and position of the break energy; 
the Band parameters are shown in the left lower part of the figure. Note, that the adopted v value would imply /3 = —2.42; however, 
relatively narrow region of the high-energy part of the spectrum and also added noise result in flatter asymptote with /3 Si —1.77. 



randomly other involved parameters from the parent normal distributions described in the previous section, we generate 
around 5,000 individual DSR spectra, add noise to them at the level of 25%, and fit each of them to the phenomenological 
Band function, see an example in Figure (9] This yields distributions of the spectral fitting parameters to be compared with 
the observed histograms of the Band parameters. 

Figure [To] displays the model histograms. Remarkably, the model -Bbrcak histogram is very similar to the observed one, 
while the a and /3 histograms display the peak values consistent with the observed ones (—1 and —2.2 respectively), although 
the widths of the model histograms are much smaller than of the observed ones. This inconsistency can be related to (i) use 
of the simplified perturbative treatment, (ii) non-optimal parameter range, or (iii) fundamental shortage of the adopted DSR 
model. We, thus, address these issues by applying full non-perturbative DSR treatment and exploring more complete range of 
the involved parameters. The remaining residuals between the model and observations will then be critically discussed within 
simplifications and limitations of the model used. 
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Figure 10. Histograms of the Band parameters a, 0, and -Ebrcak obtained within the perturbativc (jitter) DSR model assuming the 
random magnetic field is weak. 
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Figure 11. Histograms of the Band parameters a, 
magnetic field is weak, A = 1. 



and Bbrcak obtained within the non-perturbative DSR model assuming the random 



3.3 Case II: Week random magnetic field; non-perturbative treatment 



As is known i Toptygin fc Fleishnianlll987l : lFleishnianll200q ) non-perturbative treatment may be required even for tlie case of 
relatively weak random field if the energy of radiating electron is large; it results in the ex uj~^'^ asymptote, see, e.g., Figure[Sl 
It is quite clear that an immediate outcome of this non-perturbative asymptote is appearance of a values around —1/2, so 
the range of a values between —1 and —1/2 will be filled. The distribution of individual a values in this range will change 
depending on the adopted A = ujq /ujst value. 

Figure [11] displays an example of the Band parameter distributions obtained within the non-perturbative treatment 
for the case A = 1, wher e the energy content of the magnetic field es (and, accordingly, e^, to keep eg = y/^, see, e.g., 
Sironi fc Spitkovskvll2009l . and references therein) is increased to es = 10~^ to keep the break energy within the required 
window of observations. We see that the peak of the a histogram shifts to the value « —0.6, related to the ex: w"^ " asymptote. 
This histogram remains very narrow like in the perturbative case, and, furthermore, its peak value does not agree with the 
observed peak value any longer. Variation of parameter A within the range corresponding to the weak field case, A > 1 does 
not improve the situation: the distribution remains narrow with the peak value between —1 and —0.5. We conclude that the 
DSR model with the weak random magnetic field, either perturbative or non-perturbative, cannot offer a consistent fit to 
the observed a histogram, while the model -Ebrcak histogram matches well the obser ved one. It is worthwhile to note here 
that a model with a weak random field was also criticized from another p erspective llKumar fc McMahorJ 120081 ) as it may 
imply an unrealistically high level of inverse Compton emission. In addition. iKirk fc Revilld (|201G ) argued that the weak- field 
case (needed to mediate the jitter-like regime of DSR) seems to be in contradiction with the required high efficiency of the 
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Figure 12. Example of Gamma distribution of the spectral index u used in the modeling. 



particle acceleration at the shocks and strong magnetic fluctuations are needed to self-consistently accelerate electrons up to 
the gamma-ray producing energies. Thus, we turn now to the case of the strong random magnetic field. 



3.4 Case III: Strong random magnetic field 

Having the weak random field model (jitter regime) rejected, we turn now to analysis of the strong random field case, A < 1. As 
has been explained in i; 13.11 new asymptotes, including (x uj~^'^^^'''^, arise in this case, which can yield broader a distribution 
if this new asymptote comes into play. Therefore, the model a distribution will depend on adopted v distribution, which is 
in fact unconstrained by the observations. However, we can take advantage of the fact that in our model the /3 distribution 
is straightforwardly determined by the u distribution. Indeed, because /? = —v — 1, we can simply derive the required v 
distribution from the observed /3 distribution. 

The observed /? histogram reaches the peak at around —2.2 and has asymmetric skew shape with a longer tail towards 
smaller values, which apparently cannot be described by a symmetric normal distribution adopted above. Thus, we have 
to adopt another reasonably simple distribution roughly matching the observed one. Specifically, we find that the Gamma 
distribution 



F{u) = {u^ voY 



bRT{R) ' 



(20) 



with R — 2, where T{R) is the Euler Gamma function, is well suited for our modeling. This distribution has a peak at 
u — vo + b; we expect to obtain correct /3 histogram for uq — 0.9 — 1 and 6 — 0.3 — 0.2, see Figure [12] 

Other than the v distribution adopted to obey Gamma distribution with the specified parameters, all other steps of the 
modeling are the same as before. We performed many runs changing the A value and also varying other involved parameters 
within the adopted limits and found that typically the a distribution is much broader than in the case of the weak field 
considered in two previous sections in a general agreement with observations; the peak value of the a histogram varies 
between —1.1 and —0.8 when A < 1. Figure [13] displays an example of the model Band spectral parameters obtained for the 
strong random field case. 

The model results are in a remarkable agreement with the observations. Indeed, the a histogram is a symmetric one, it 
displays a peak at the right place, a = —1, and its bandwidth is comparable to that of the observed histogram. The /? histogram 
almost repeats the observed one, displaying the correct asymmetric shape and the peak at the right place, /? = —2.2. The 
inbreak histogram agrees with the observed one rather well: it has correct shape and bandwidth, although the peak value is less 
than the observed value by the factor around 2. This discrepancy is, however, inessential: as we have seen in ij 13.31 the position 
of the break energy can easily be adjusted by a small change of the magnetic energy content es and corresponding change of 
ee. Moreover, considering broader range of ^e variation (recall, we adopted a constant value of ^e = 1 in our modeling) can 
also easily change the characteristic break energy by a factor of 2 or more. Thus, we can conclude that the simplified DSR 
model adopted in this section is intrinsically capable of reproducing the Band parameter distributions compatible with the 
observed ones. 
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Figure 13. Example of the model Band parameter (a, /3, and Ebreak) distributions obtained within the DSR model with strong random 
magnetic field. 
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Figure 14. Cross-correlation of tlie model Bend spectral parameters a and /3. 
Figure 15. Cross-correlation of the model Bend spectral index a and Ep^^k = (a -I- 2)Eo. 
Figure 16. Cross-correlation of the model Bend spectral index /3 and -Epeafc = {ct + 2)Eo. 



3.5 Cross-correlations between the model Band parameters 

In addition to the Band parameter distributions themselves, it is worthwhile to address a question if our model reproduces 
the cross-correlations between the Band parameters correctly. This task can easily be solved using the variety of the model 
spectra produced in each model run. 

Figures [1411161 displav these cross-correlations to be compared with Figure 31 from iKaneko et al.l ( 20061 ). Like in the 
observation, the spectral indices a and /? are not hi ghly correlated, alth ough in the model plot the region of —0.5 < a < 
is underpopulated compared with the observed plot ( Kaneko et al.l 20061 ) . Two other plots are in remarkable agreement with 
the observed cross-correlation plots, presented in iKaneko et al.l ( 20061 ). We conclude that the developed model is naturally 
capable of reproducing the cross-correlation plots in addition to the histograms themselves, which is a remarkable success of 
the non-perturbative DSR model in the presence of strong random magnetic field. 



4 DISCUSSION 

Our modeling shows that we can get an overall agreement between the model and observed histograms of the Band GRB 
spectral parameters within the non-perturbative DSR model with strong random magnetic field when A = uon/uost ~ 0.015. 
We note that in our model the A parameter is not derived from a microscopic treatment of the shock interactions, rather it 
is a free model parameter adjusted for the model histogram to resemble the observed ones. Let us consider if the obtained 
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Figure 17. Histogram of the ujQ/ujpB ratio derived from tfie DSR model with the strong random field with A = 0.015. 



A value has any sense versus current models of magnetic field generation at the shock fronts. To do so we recall that cjq is 
defined by the correlation length Lmax of the random magnetic field, cjq ~ 27rc/Z/max. If the magnetic field is produced by a 
two-stream (filamen tation or Weibel) instabilit y, the correlation length Lmax is expected to be about ten plasma skin scales. 



10c/a;po (e.g., ISironi fc Spitkovskvll2009l ). Since the frequencies ujpc and ljq are independent input parameters in the 



L 

modeling, we can check if the correlation length is indeed of the order of the skin scale, by inspecting the actual distribution 

of the ujo/ujpc ratio as it appears in the best-fit model case. 

Figure [17] displays a histogram of the cjo/i^pc ratio obtained for the model with A — 0.015. The distribution has a 
symmetric bell shape with the peak about 0.5, thus, Lmax ~ 4:TtIsc and the required random field correlation length is indeed 
of the order of ten plasma skin scales, which agrees with the idea of the random field generation by a two-steam instability in 
the internal shock interactions. 

We must note that although the model histograms and the cross-correlation plots look similar to the observed ones, there 
are some residuals between them. For example, the model a histogram is somewhat narrower than the observed one with a 
deficit of values —2 < a < —1.5 and —0.5 < a < 0. The first interval can possibly be filled if one considers a nommiformity of 
the emission source, which can naturally broaden the emission spectrum leading eventually to smaller a values. The second 
interval requires some additional physics, not included in our simplified modeling, to be taken into account, for example, 
random electric fields, turbulence anisotropy, or specific source geometry/viewing angle combination. 

On the other hand, given the number of simplifying assumptions adopted for the modeling, we can conclude that the 
obtained agreement between the model and observations is remarkably good. Let us briefiy remind and discuss those simpli- 
fications. 

(i) We considered random magnetic fields only. I t is known, however, that diff usive radiation in electrostatic Langmuir 
waves results in a harder spectrum, e.g. I^, oc uj^ , ( Fleishman fc Toptvginll2007al la). thus, the presence of these Langmuir 
waves could compensate the deficit of the —0.5 < a < values. 

(ii) A monoenergetic spectrum of accelerated electrons was adopted. Although typically the distribution of /3 parameter is 
ascribed to a parent distribution of the spectral index p of the electron distribution over energy, N{E) ex E^ , our modeling 
shows that the right distribution of the /? index can easily be obtained even for a monoenergetic electron spectrum N{E) oc 
5{E — £*). In fact, for power-law energy distributions almost the same results hold for p > 2i/ + 1. However, if a power-law 
range with p < 21/ + 1 is present in the electron energy spectrum, this is not in a contradiction with the model, although this 
does add more flexibility to formation of the /3 distribution, which can further broaden it towards even better similarity to 
the observed histogram. 

(iii) A number of the source parameters (e.g., z and ^e) were adopted to be the same for all the sources; in fact, adopting 
more realistic distributions can broaden the obtained distributions of the Band spectral parameters and also affect the best-fit 
parameter space, so the accuracy of the obtained source parameters is at best to an order of magnitude. 

(iv) Random magnetic fields accounted by the model were adopted to be statistically uniform and isotropic. Inclusion of 
the turbulence anisotropy can affect the radiation spectra and so modify the Band parameter distributions. 

(v) Having adopted both turbulence and electrons are isotropically distributed, we did not consider any source geome- 
try/viewing angle effect. For the case of anisotropic distributions such effects can also come into play. 
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(vi) And finally, we did not explicitly consider the source evolution, although both particle distribution and magnetic field 
can evolve in time resulting in a GRB spectral evolution. 

The latter three simplifications have in fact been addressed (e.g.. lWorkman et al.ll2008l . and references therin) in a number 



of studies, however, a purely perturbative (jitter) weak-field regime of the DSR was adopted to calculate the radiation spectra. 



This weak-field regime wa s shown to rise severe problems in application to the GRB prompt emission (JKumar fc McMahon 



200a : lKirk fc Revillell20ld ). In addition, as has been demonstrated above in this paper, this jitter regime of the DSR is incon- 



sistent with the observed Band parameter distributions, while the strong-field regime is in fact needed. Since the strong-field 
regime requires a much more sophisticated fully non-perurbative treatment, those previous studies cannot be straightforwardly 
applied to this case, and so the required generalization must be specifically performed from scratch within the non-perturbative 
treatment iJToDtvgin fc Fleishmanlll987l : lFieishman|[20oi : [Fleishman fc BietenholJboOTl '). 



Although all these effects are potentially important and must eventually be taken into account in building a more 
comprehensive model, we conclude that even the simplified DSR model considered here is naturally capable of reproducing 
main characteristic properties of the Band parameter distributions and cross-correlations between them. The ranges of the 
parameters needed for the model to most closely reproduce the observed histograms agree well with standard fireball model 
parameters. 
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